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ABSTRACT 

We present high resolution (R — 75,000-100,000) mid-infrared spectra of the high-mass embedded 
young star IRS 1 in the NGC 7538 star-forming region. Absorption Hnes from many rotational states 
of C2H2, 13C12CH2, CH3, CH4, NH3, HCN, HNCO, and CS are seen. The gas temperature, column 
density, covering factor, line width, and Doppler shift for each molecule are derived. All molecules 
were fit with two velocity components between -54 and -63 km s~^. We find high column densities 
(~ 10^^ cm~2) for all the observed molecules compared to values previously reported and present 
new results for CH3 and HNCO. Several physical and chemical models are considered. The favored 
model involves a nearly edge-on disk around a massive star. Radiation from dust in the inner disk 
passes through the disk atmosphere, where large molecular column densities can produce the observed 
absorption line spectrum. 

Subject headings: Infrared: ISM - ISM: molecules - HII regions: NGC 7538 IRS 1 



1. INTRODUCTION 

The NGC 7538 (S158) molecular clo ud, located in the 
Perse us arm at a distance of 2.8 kpc (jUngerechts et all 
|2000( ). contains many protostellar objects, making it a 
good candidate to study different stages of star for- 
mation. Furthermore, there seems to be a progression 
in evolution from northwest (the visible H II region, 
NGC 7538, and IRS 2, an only moderately obscured com- 
pact HII region) to southeast (IRS 9, a deeply embedded 
core) (jElmegreen fc Ladal [19771 : iCampbell fc ThomDSonI 
Il984[ ). The near- and mid- infrared sources IRS 1, 2, and 
3 all lie near the northwest end of the star-forming region. 
Of these, IRS 1 is the brightest r nid-infrared and radio 
continuum source (iMartinI [1971 |W vnn- Williams et al.l 
119741 : rWillnerj[T97^ . Its exciting star has a luminosity of 
L > 8 X 10** Lq, which corresponds to a ZAMS spectral 
type earlier t han 07.5, and an ionizing photon fl ux of 
> s-i (|Werner et al.l[l979l : iLugo et al.l[200l . 

IRS 1 is embedded in a dense molecular cloud. Mil- 
limet er continuum emission is seen from t he sur rounding 
dust ([Scoville et al.lll986l ). iPratap et aP (|1989D mapped 
the HCO+ and HCN emission from within ~1' of IRS 1 
with 3" resolution. They found evidence f or a shell-like 
struct ure on scales of 10-20" around IRS 1. IWilson et al.l 
(Il983f) observed NH3 absorption toward IRS 1, and 
iHenkel. Wilson. &: Johnstonl ()1984D mapped the absorp- 
tion with the VLA. The absorption traces warm gas 
with a high column density of NH3 and is centered near 
—60 km s^^. The emission, on the other hand, is cen- 
tered near —56.5 km s~ ^. The systemic velocity for this 
source is ~ 57 km s~ (Ivan der Tak et a"ni2OOO0 . More 
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recently, IZheng et al.l ()2001[ ) mapped NH3 emission from 
a more extended region in NGC 7538. This emission 
probably probes the outer envelope as indicated by the 
colder temperature and the smaller column density. Us- 
ing both single-dish and interferometer observati ons of 
various molecular tracers . Ivan der Tak et all (|2000D stud- 
ied the density and temperature structure of both the 
cold outer envelope and the warm inner material (240- 
72000 AU). The inner warm region is characterized by 
temperatures of a few hundred K. 

Outflows are often present to ward protostars, and 
NGC 7538 IRS 1 is no exception. iCampbelll (fl98l ob- 
served IRS 1 with the VLA at 5 and 15 GHz. She found a 
pair of very compact lobes of continuum radiation, sep- 
arated in declination by 0.2", with emission extending 
out to ±2". Her preferred model involves a bipolar ion- 
ized outflow from a late O star, coUimated by a core of 
dense gas extending from <65 AU to >25, 000 AU. Fur- 
ther evidence of outflows was found when Gaum e et al.l 
([1993) observed IRS 1 with the VLA in the H66a line 
and the 22 GHz continuum. Their continuum image is 
similar to that of Campbell, but shows additional struc- 
ture. They found broad recombination line emission in 
the two lobes, with a minimum line-to-continuum ratio, 
which they attribute to high electron density, between 
the lobes. They propose a model involving a high veloc- 
ity stellar wind interacting with photo-evaporating knots 
of neutral gas. 

More recently, iLugo et "aD (J2004') have modeled the ra- 
dio continuum observations of IRS 1 as due to a wind 
produced by photo-evaporation of a circumstellar disk 
with a radius of 500 AU exposed to UV radiation from 
the central O star. M aser lin e emission is seen from H2O 
(iKamcva et al.' 1990^, OH (Dickel et al.1 [l982il. H 9CO 
Ikoffman et al. 2003} CH3 OH (Me nten et al.lll986D. and 
NH-i (iMadden et all lijSl . iMinier. Booth, fc Conwavl 
()1998( ) used VLBI observations of the methanol masers 
to infer the presence of a nearly edge-on rotating disk 
traced by one of the maser clusters toward IRS 1 . The 
central velocity of this cluster, which is the brightest 
cluster near IRS 1 in both CH3OH and NH3 masers, 
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is -56.2 km s^^. A second cluster of masers, 0.25" to 
the south, is at -61.0 km s~^. The disk masers he along 
a line centered near the gap between the lobes of the 
free-free continuum radiation, although they are aligned 
along a position angle ^^30° from the symmetry plane. 
iPestalozzi et al.l (|2004l ) made a detailed model of maser 
emission f rom a disk ar ound IRS 1, supporting the sug- 
gestion of iMinier et al.l They concluded that the maser 
emission comes from disk radii of ~ 290 —J50 AU . How- 
ever, more recently, iDe Buizer fc Minieil ()2005D show 
that the CH3O H masers may trace k nots in the outflow 
instead. Also, iDe Buizer fc Minierl identify a disk-like 
structure perpendicular to the outflow direction in mid- 
infrared continuum. The extent of the infrared emission 
is -450 AU. 

Infrared absorption by dust grains and their icy man - 
tles has also been observed toward IRS 1. iWillneii (fT976h 
and lWillner et al.l (|1982f) observed 9.7 fiui silicate absorp- 
tion and absorption by various ices, including H2O. The 
ice bands toward IRS 1 are relatively weak, especially in 
comparison with NGC 7538 IRS 9, which must l i e in a 
colder region of the NGC 7538 cloud. iGibb et all (l200i ) 
observed the complete 2.5-20 fim region with the Short- 
Wavelength Spectrometer (SWS) on the Infrared Space 
Observatory, ISO. II2O, CO, and CO2 ices are clearly 
present, as well as other less securely identified features. 
Toward NGC 7538 IRS 9, -1' southeast of IRS 1, a much 
richer ice spectrum is seen, with CH3OII, XCN (probably 
OCN~), and CII4 also confidently identified, and OCS, 
H2CO, HCOOH, and NH3 hkcly present. 

Gas-phase absorption in the infr ared has been observed 
toward IRS 1 by several groups. [Mitchell et al.l ()1990D 
observed absorption by ^^CO and ^^CO in their 5 /im 
V = 1 — bands. They did not spectrally resolve the 
lines, and they used a curve of growth analysis assum- 
ing pure absorption by Gaussian lines with FWHM of 
8 km s""'^ to derive the gas temperature and column den- 
sity from the less-saturated ^^CO lines. They found two 
absorbing components: cold gas at 25 K, and warm gas 
at 176 K. A gas density of nn > 10^ cm~'^ was required 
to maintain the population of t he high- J levels in the 
warm gas. Using the ISO-SWS, iLahuis fc van Dishoeckl 
(|2000D detecte d warmer gas than the CO observed by 
iMitchell et al.l . They observed the 1/5 band of C2II2 at 
13.7 urn and the 1^2 band of HCN at 14.0 /xm. With 
the SWS resolution of —1800, they were not able to re- 
solve individual line shapes, and were most sensitive to 
the blended lines of the Q branches, but they were able 
to derive temperatures and column densities from the 
shapes of the Q branches and the depths of the v = 1 — 
fundamental and u = 2 — 1 hot band Q branches. They 
derived T 800 K and TV = 0.8 x lO^*^ cm'^ for C2H2, 
and T ^ 600 K and N = 1.0 x 10^^ cm'^ for HCN. 
iBoonman efall (|2003f l derived T = 500 K for C2H2 from 
an update d reduction of the spectra . |B oonman_ et al.l 
([2003) and IBoonman fc van Dishoeckl ((2003) studied gas 
absorption from CO2 and H2O using ISO data. They 
found enhanced abundances toward the inner warm ma- 
terial compared to the cold envelope indicating that grain 
mantles are sublimating and enriching the gas-phase 
chemistry close to the pr otostar. Gas-phase CH4 was also 
detected toward IRS 9 (|Lacv et al.l 119911 : iBoogert et al.l 
HqM) but not IRS 1. 

The wealth of observational data available suggest the 



following physical scenario for NGC 7538 IRS 1. The 
massive young star is surroun ded by a massive cold en - 
velope on scales - 72000 AU (|van der Tak et al.ll2000l ). 
At scales below 1000 AU, the temperature increases to 
> 100 K. Inside this radiu s, there is evidence for a disk 
( De Buizer fc MinieH l2005f ) and small scale knotty out- 
flows ( Gaume et al.lll995i r both of which can be affected 
by the UV and X-ray rad iation from the protostar (e.g., 
iStauber et a"ni2003 . l2005h . 

In this paper, we present high resolution mid-infrared 
spectroscopy of NGC 7538 IRS 1 showing a rich ab- 
sorption spectrum containing lines from seven molecules: 
C2H2, HCN, CH4, CH3, NH3, HNCO, and CS. Mid- 
infrared absorption is much more sensitive to the in- 
ner warm gas than radio observations. Previous infrared 
studies suffered from limited spectral resolution. Section 
m describes the observations. We then present the model 
used to derive the column densities and temperatures for 
the various species in section [3] and describe the results 
for each of the molecules in section U) Subsequently, in 
section jSJ we discuss the possible scenario in which we 
are probing chemistry in a circumstellar disk. In section 
|6l we provide some concluding remarks. 

2. OBSERVATIONS 

2.1. Observations and Data Reduction 

The observations were made with TEXES, the Texa s 
Echelon Cross Echelle Spectrograph ()Lacv et al.ll2002f) . 
on the NASA Infrared Telescope Facility (IRTF) in 2001 
June and November, 2002 September and December, 
and 2005 December. TEXES is a high-resolution cross- 
dispersed spectrograph operating at mid-infrared wave- 
lengths between 5 and 25 fim. It achieves a spectral 
resolution of i? = 75,000-100,000 (Aw = 3 - 4 km s^^) 
shortward of 14 /im with a slit width of ~1.5" . (Of 
the observations presented here, only the spectra near 
11 fim, which used a 1" slit, achieved the higher resolu- 
tion.) The spatial resolution along the North-South ori- 
ented slit is ~1", slightly larger than the diffraction limit 
of the IRTF. The spectral and spatial sampling by the 
256 X 256 pixel Si: As detector array are 1.0 km s~^ and 
0.35", respectively. At each spectral setting, 5-10 orders 
of the high resolution echelon grating are recorded, giving 
a spectral coverage of AAtot /A ~ 0.5%. At wavelengths 
shortward of 11 /im full order widths are observed, giv- 
ing continuous spectral coverage, but longward of 11 /im 
the order width exceeds the array width, leaving gaps 
in the observed spectra. The location of the gaps can be 
shifted in order to optimize line coverage. However, since 
the gaps are small and, at 13 /im, there are many telluric 
lines, it would be impractical to repeat observations with 
small shifts in order to get complete spectral coverage for 
each setting. 

Much of the data were taken in high water vapor con- 
ditions, but the observing conditions did not seriously 
affect the spectra except for increased noise, especially 
on atmospheric lines. Observations of NGC 7538 IRS 1 
were interspersed with observations of a bright asteroid 
(usually Ceres (> 500 Jy) , but also Hygeia (~ 150 Jy), 
which served as a comparison source for removal of tel- 
luric absorption features. The telescope was nodded by 
3-5", moving the source along the spectrograph entrance 
slit to allow subtraction of background emission. At the 
beginning of each set of nodded observations an ambi- 
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ent temperature blackbody was observed for flat-fielding. 
The first nods were used to peak up on the source by 
maximizing the throughput of the sht. The nod pairs 
taken during peak up were not included in the final sum. 

Data reduction followe d the standard TEXES pro- 
cedure (jLacv et al.l l2002[ ) of first subtracting the two 
nod position spectra, flat- fielding with the blackbody- 
sky difference echellogram ()Lacv et al.lll989D . and inter- 
polating over spikes and bad pixels. Optimally-weighted 
point-source spectra were then extracted from the echel- 
lograms. The spectra were linearized in In(^) based on 
the known optical distortions, and absolute wavenumber 
calibration was obtained from telluric emission features 
and corrected for the Earth's motion relative to LSR. The 
resulting wavenumber scale is correct to within 1 km s~^. 
The same procedure was used for the asteroid spectra, 
and then the IRS 1 spectra were divided by the asteroid 
spectra with correction for differences in airmass. With 
this procedure, telluric absorption fines as deep as 90% 
can be removed, although of course the noise increases 
on lines. However, broad response variations often re- 
mained, which were removed by fitting the continuum in 
each echelon order to a quartic polynomial and dividing. 
This procedure also sets the continuum in each order to 
one. 

2.2. Line Identification and New Infrared Detections 

Ten spectral settings were observed in the 728- 
820 cm-i (13.7-12.2 ^m) region, two in the 860-930 
cm~^ (11.5-10.7 ixm) region (see Figure [1]), and three 
in the 1240-1312 cm~^ (8.0-7.6 /xm) region (see Figure 
[2]). A total of 110 echelon orders, or ~70 cm~^, were 
observed. The original observations were meant to study 
the line profiles of C2II2 and HCN, which had previously 
been detected. In addition to the prominent fundamen- 
tal band lines of C2H2 and HCN, many C2H2 hot band 
lines were detected in our spectra including lines from the 
+ f 5 - and the 21/5 - 2/5 bands. After the hot bands 
were identified, unidentified lines still remained, some of 
which were double dipped. The double lines turned out 
to be CH3 lines which are split due to spin-rotation inter- 
actions (see Appendix). This is the first detection of CH3 
toward dense gas and the first ground-based detection of 
CH3. Previously, CH3 had only been observed with ISO 
toward the Galactic center (|Feuchtgruber et al.ll2000l ) . In 
efforts to detect ethane, C2H6, at 810-820 cm~^, many 
lines were detected, though none corresponded to C2H6. 
Some were due to NH3, but many lines remained uniden- 
tified (see Figures [3] and H]). Based on their separation 
we were able to identify the lines as due to HNCO. In- 
terstellar HNCO is a well kn own h ot cor e molecule at 
radio wavelengths ijZinchenko et all l2000f ) , but it had 
never been seen at infrared wavelengths. In total, five 
molecules were observed in the 728-820 cm~^ region: 
C2H2 (including ^^C^'^CYL2), HCN, HNCO, CH3, and 
NH3. At higher frequencies, CH4, C2H2 and CS were ob- 
served in the 1240-1320 cm~^ region. CH4 observations 
from the ground are possible for this source because of 
the large Doppler shift with respect to the telluric lines. 

Table [T] gives a summary of the detected lines. At 
least 5 lines were observed for each molecule, with some 
molecules like C2H2 (including isotopic lines and hot 
bands), and NH3 having more than 40 lines detected. 
However, the number of HNCO lines dominates over 



TABLE 1 

Summary of Observed Lines 





Total Niimhpr of 
DstrGctsd Liiiss 


Band 


C2H2 


25 




C2H2'' 


2 


1^4 + i^S 


C2H2'' 


20 


Vi + I'S-fi 


i3ci2cH2'= 


7 


^5 


HCN 


10 


V2 


NH3 


64 


V2 


HNCO 


125 


1/4 


CH3 


11 


Vl 


CH4 


12 


U2/VH dyad 


CS 


6 


V2 



Note. — a. These lines were observed at 8 

b. Lines in Q-brancli at 731 cm~^. See top 
panel of Figure [T] 

c. Only unblended lines are reported here. 



TABLE 2 
Observed Lines 



Molecule 


Line 


Wavenumber 


Ej 


Band 






cm~^ 


(K) 




C2H2 


Q(5) 


729.289 


35.30 




C2H2 


Q(6) 


729.343 


49.42 


1^5 


C2H2 


Q(7) 


729.406 


65.89 


1^5 


C2H2 


Q(8) 


729.477 


84.71 


1^5 


C2H2 


Q(9) 


729.558 


105.89 


1^5 


C2H2 


Q(10) 


729.648 


129.41 


1^5 


C2H2 


Q(ii) 


729.747 


155.29 


1^5 


C2H2 


Q(14) 


730.096 


247.02 


1^5 


C2H2 


Q(15) 


730.230 


282.30 


1^5 


C2H2 


Q(16) 


730.373 


319.93 


1^5 


C2H2 


Q(18) 


730.686 


402.22 


1^5 


C2H2 


Q(19) 


730.855 


446.89 


1^5 


C2H2 


Q(20) 


731.034 


493.90 


1^5 


C2H2 


Q(26) 


732.287 


825.20 


1^5 


C2H2 


Q(27) 


732.526 


888.62 




Note. - 


— Only lines that were 


individually de- 



tected at the 3-cr level are included in this list. Full 
line list is available online only. 

lines from the other molecules, with over 100 lines de- 
tected. Table[2](full version is available online only) lists 
all the individual lines detected in our spectra. A se- 
lection of spectra including all observed molecules are 
shown in Figures [T] and [H along with fits described 
in For most molecules the line strengths were 

taken from the GEISA03 dat abase^. The CH3 band 
streng th was obtained from IWormhoudt fc McCurdvl 
(11989') and the HNCO oscil lator strength was obtained 
from Lowenth al et al.l (12003). The CS oscilla tor strength 
was taken from iBotschwina fc Sebaldl (|1985f ). 

3. SPECTRUM FITTING 

To derive physical parameters of the absorbing gas, 
model spectra were fitted to the observations. For this 
purpose, the data were divided into two groups. All data 
at wavenumbers between 700 and 1000 cm~^ (10.7-13.7 
/im) were fitted simultaneously, and data at 1240-1280 
cm~^ (near 8 /im) were fitted s eparat ely. The Mar- 
quardt fitting procedure (|Bevington fc RobinsonI [20031) 
was used, which minimizes x2, the summed squared de- 

http:/ /ara. lmd.polytechnique.fr 
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729.5 



730.0 



730.5 



731.0 



X 



CD 
N 




II II CH, 



762.0 



762.5 



763.0 





807.5 




1.1 
1.0 

0.9 

0.8 

0.7 
866.5 













NH3 1 


1 1 II \ 



867.0 867.5 868.0 

Wavenumber (cm~^) 



868.5 



Fig. 1. — Selected orders (thick black) between 729 and 870 cm ^ are shown with the best fit (red). Positions of lines considered in the 
fit are given: C2H2 (light blue), l3cl2CH2 (blue), hot bands of C2H2 (dark blue), HCN (green), CH3 (orange), HNCO (red), and NH3 
(purple). The thin black line indicates the atmospheric transmission. The top panel has a larger scale for the normalized fiux (y-axis) to 
show the lower atmospheric transmission. Note that there are gaps between orders. 
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1245.5 1246.0 1246.5 1247.0 1247.5 




1248.0 1248.5 1249.0 1249.5 1250.0 1250.5 



o 




1278.0 1278.5 1279.0 1279.5 1280.0 




1282.0 1282.5 1283.0 1283.5 1284.0 

Wavenumber (cm~^) 



Fig. 2. — Selected orders (thick black) between 1245 and 1285 cm~^ are shown with the best fit (red). Position of lines considered in 
the fit are given: CH4 (blue) and CS (red). The thin black line indicates the atmospheric transmission. In this region there are no gaps 
between orders. Instead the orders overlap where the black line is darker. 
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: HNCo 




I I 



805 810 815 820 

Wavenurrber [cm"'] 

Fig. 3. — This figures shows a portion of spectrum containing 
some NH3 lines as well as wealth of HNCO lines. These spectral 
settings led to the identification of HNCO. Figure |4] shows a small 
portion of this spectrum. 



Warm Continuum 
Source 




Absorbing 
components 



> 



Observer 



Fig. 5. — This figure shows two absorbing components between 
the continuum source and the observer. Each component covers 
10% of the continuum source. The overlap of the two absorbing 
regions is given by the product of their covering factor. In this case, 
only 1% of the continuum source is covered by both components. 




I HNCO R{10) 
: I1U-IU 
' 7 I HNCO R(10) 
- ML-1L 

NHjl 
P(8,6)sl 



I HNCO R(1 1) HNCO R(12) 

I1U-IU 1U-1U 

|HNC0R(11) HNCO R( 12)1 NHjl 

1 1L-1L 1L-1L I P(8,1)sl 

NH3I NH3I NH3I 

P(8,4)s I P(8,3)s I P(8,2)s I 



o.a [ L 

806.0 



Fig. 4. — This figure shows the profiles of several HNCO and 
NH3 lines in the spectral region shown in Figure|3] The two HNCO 
sub-bands shown are R-branch transitions where the K = 1 state 
is split into upper (U) and lower (L) levels. Because of the close 
separation the sub-bands, the lines can overlap for a range of J. 
See Appendix A for more details on transition rules for HNCO. 



viation normalized by the squared noise, between the 
data and the model. 

Fitting was first attempted with a single component 
for each molecule, but the residuals suggested that a 
better fit could be obtained with multiple velocity com- 
ponents. The model that we used assumed that each 
molecule is found in one or two absorbing components. 
Each component is described by its Doppler shift, Vlsr, 
its Gaussian line width (1/e half width = (2kT/m)i/2), 
b, the molecule's column density, N, the rotational tem- 
perature of the gas, T, and the covering factor, C, (the 
fraction of the background continuum source covered by 
the component). The fitted column density is the av- 
erage over the partially covered source, so the column 
density in the covered portion would be N/C. A cov- 
ering factor of, e.g., 10% means that lines saturate at 
90% of the continuum flux and the column density in 
the absorbing component is ten times the average over 
the continuum source. The effect modeled in this way 
could also result from veiling (continuum emission from 



foreground or surrounding material) or from re-emission 
in the lines by the absorbing molecules. The components 
were assumed to overlap by the products of their cover- 
ing factors (see Fig. [5|) . The following equation was used 
to calculate the observed transmission: 

lobs = /o(l - Ci(l - e-"i('^)))(l - C2{1 - e^^^f"^))), (1) 

where Ci and C2 are the covering factors for the two 
components. This results in a greater absorption in sat- 
urated lines than would result if the optical depths were 
added first and then the transmission spectrum was cal- 
culated from the optical depth spectrum. This approach 
was chosen because it gave a better fit to the data than 
adding the optical depths first and assuming the same 
covering factor, but it assumes that different components 
absorb along different lines of sight to a partially covered 
continuum source, which may not be the case. 

A frequency correction was allowed for each spectral 
setting observed, to correct for errors in wavelength cali- 
bration. With only one exception, the correction was less 
than a 1 km s^^ Doppler shift. In addition, a broadening 
of the instrumental resolution was permitted in the fit- 
ting to allow for imperfect internal instrument focusing, 
or possibly a 'macroturbulent' broadening in the absorb- 
ing gas. The fitting program chose a resolution a factor 
of 1.2 greater than that derived from gas cell data. A 
constant frequency resolution (as opposed to a constant 
Doppler resolution) was used in each of the two fitted re- 
gions. This represents well the improving resolving power 
toward shorter wavelengths in the 10.7-13.7 fim region. 
In addition to these parameters, the continuum, slope, 
and curvature of each echelon order were varied, to allow 
correction of the baseline fitting done during data re- 
duction. The continuum fitting and frequency correction 
required about 300 free, but rather easily determined, 
parameters. Fewer parameters are needed to determine 
the physical conditions such as temperature, column den- 
sity, line width and covering factor: 70 parameters for the 
10.7-13.7 /im spectra and 21 for the 8 /im region. For the 
10-13 /im region, we have over 3000 points to constrain 
the 70 parameters of interest if we characterize the con- 
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straining points by the number of lines times the number 
of pixels for each line. 

For C2H2 and HNCO, lines of several bands were ob- 
served. For C2H2, in addition to lines of the Q and R 
branches of the fundamental, lines of several Vi + - 
vc, and 2i^5 - z/5 bands (involving absorption from the ex- 
cited and i/5 vibration states) were observed. The vi- 
brational temperatures describing the populations of the 
and 1^5 states were included as free parameters in the 
fit. R-branch lines of the vc, band of ^•^C^^CH2 were also 
observed, and the ^"^0^^0112/02112 abundance ratio was 
allowed to vary, although it was kept the same in all the 
O2II2 components. Two O2H2 lines of the z/4 -I- 2/5 band 
at 7.6 /im were observed as well. For HNOO, lines from 
the band were observed. HNCO is an only slightly 
non-linear molecule, and its spectrum consists of a se- 
ries of sub-bands resembling those of a linear molecule, 
like O2II2. Its dipole moment oscillates diagonally, giv- 
ing it a- type sub-bands, with no change in the angu- 
lar momentum about its long axis and only P and R 
branches, and b-type sub-bands, in which the angular 
momentum about its long axis changes, and P, Q, and R 
branches are seen. See Appendix for further discussion 
of the HNOO spectrum. Since we did not have labora- 
tory data regarding the relative strengths of the different 
sub-bands, we considered them to be caused by different 
species and summed the abundances to obtain the HNOO 
abundance. The HON lines observed were from the V2 
bending mode. We included HON 2v2 - V2 lines in the 
fit. Although no hot HON lines were obvious in the spec- 
tra, their inclusion lowered significantly (Ax^ — 13) 
and changed the best-fitting A'hcn noticeably. 

Although it might not be apparent from the spectra, 
given the shallowness of the lines, many of the detected 
lines are at least moderately saturated, and some are very 
saturated. The lines only appear weak because of the 
small covering factors (or emission filling in the absorp- 
tion lines) and to a lesser extent the small line widths. 
This conclusion is very robust; the relative depths of lines 
of different opacity requires substantial saturation. This 
is most apparent in the near equality of the depths of 
the ortho and para lines in the O2H2 Q branch, near 730 
cm~^ (see top panel of Figure[T]). HON in the -60 km s~^ 
component and NH3 in the -56.5 km s~^ component are 
so saturated that only lower limits could be placed on 
their abundances until their isotopomers were included 
in the fit to provide upper limit constraints. By including 
isotopomers and intrinsically weak lines, a wide enough 
range of line optical depths was observed to allow mean- 
ingful constraints to be placed on the molecular abun- 
dances in spite of the saturation. 

The observed lines of OH3 were in the R branch of the 
V2 out-of-plane bending mode of this planar molecule. 
Each of these lines is doubled by the interaction between 
the electron spin and the molecular rotation (see Ap- 
pendix). The OH3 V2 band center is at 16.5 /im, whereas 
all observed lines were shortward of 13.5 /im, with J > 
7, making us rather insensitive to cold OH3. Lines of the 
P and Q branches of the NH3 1/2 band were observed. 
There was no evidence of ^^NHa lines, but they were 
included in the fit to constrain the NH3 abundance, as 
was done with H^'^ON. In the 8 /im region, which was 
fitted separately, lines of the 1^2 / 1^4 dyad of OH4 were ob- 
served. Because of the strong absorption by low- J lines 



of OH4 in the Earth's atmosphere, most of the observed 
lines originate at J > 4. We were able to observe the 
R(0) line due to the favorable Doppler shift of the source 
(~ —57 km s^"'^). The R(0) line helps constrain the tem- 
perature derived from the high J lines. We also detected 
and fitted absorption by OS. 

The derived component parameters are given in Table 
[3l Uncertainties for the parameters, given in parentheses, 
are three times the square roots of the diagonal elements 
of the error matrix, which are uncertainties allowing all 
other parameters to vary. The noise in the fitted spectra 
is assumed to be purely photon statistical noise. With 
this noise estimate, the reduced ~ 1-6, indicating that 
non-statistical noise sources contribute moderately, or 
perhaps that a different model for the line profiles is jus- 
tified. Including non-statistical noise, the uncertainties 
given might reasonably be taken to be 95% confidence in- 
tervals. The larger fractional errors are asymmetric, and 
would be more symmetric in the log of the parameters. 
For column densities, this is a result of the non-linearity 
of the curve of growth for moderately saturated lines. 
Oolumn densities with errors greater than their values are 
not consistent with zero. For unsaturated lines, the de- 
rived covering factors are very uncertain because equiv- 
alent widths depend only on the product of the covering 
factor and the column density. For example, the column 
densities of HNOO and OH3 have smaller fractional un- 
certainties than those of O2H2 and HON, whereas the 
covering factors have smaller uncertainties for O2H2 and 
HON. All of the fitted molecules are detected with very 
high confidence. 

4. RESULTS 

For most molecules, we found two components 
with Doppler shifts of Vlsr — —56.4 ±0.3 km s~^ and 
-59.7 ± 0.3 km s"!, with widths of ^ 1 - 2 km s'^ 
(see Table [SI). Most of the component temperatures were 
in the range 200-400 K and covering factors were typ- 
ically 0.04-0.08, with O2H2 requiring somewhat larger 
values of 0.06 and 0.18. HNOO has somewhat broader 
component line widths but other parameters are simi- 
lar. OH3 has the broadest component line widths, with 
a greater velocity spacing between the two components 
and the hottest component. Both OS and OH3 are fitted 
with much larger covering factors (consistent with 100%) 
than were the other molecules. However, these covering 
factors are quite uncertain due to the weakness of the ob- 
served lines, and the large values may not be significant. 
The model spectra resulting from the best fit parameters 
can be seen in Figures [1] and [2l 

In addition to agreeing among the different molecules 
(except OH3), the Doppler shifts agree with the two most 
commonly observed values for molecules seen at radio 
wavelengths. This gives us some confidence that the two 
velocity components, at -56.5 and -60 km s~^, are real, 
although whether they are two separate components, or 
just a way of describing a non-Gaussian velocity distri- 
bution is difficult to determine from the data. Gener- 
ally, the temperatures of the two components for a given 
molecule are not very different, although a fit constrain- 
ing them to be equal was significantly worse than the fit 
allowing all temperatures to vary. Line widths, column 
densities, and covering factors vary more. The fit was 
very poor when they were constrained to be equal for all 
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TABLE 3 

Results from the model fit to the spectra 



Component 1 



Component 2 



Molecule 


14SR 


b 




N{X) 


C 


Vlsr 


b 




N{X) 


C 




(km s~-^) 


(km s-i) 


(K) 


(10^6 cm-2) 


(xlOO) 


(km s"-"^) 


(km s"-*-) 


(K) 


(IOI6 cm-2) 


(xlOO) 


C2H2 


-55.7 (0.3) 


0.6 (0.1) 


225 (20) 


3.0 (0.6) 


5.9 (0.6) 


-59.4 (0.3) 


0.6 (0.1) 


191 (10) 


2.8 (0.2) 


23.7 (1.0) 


HCN 


-56.3 (0.3) 


0.8 (0.1) 


256 (30) 


5.6 (2.8) 


8.4 (0.9) 


-60.0 (0.3) 


1.0 (0.4) 


456 (127) 


1.3 (1.2) 


68.4 (49.7) 


HNCO 


-57.2 (0.3) 


1.8 (0.2) 


319 (27) 


0.4 (0.3) 


3.8 (0.3) 


-60.2 (0.2) 


0.7 (0.4) 


171 (79) 


0.1 (0.3) 


97.0 (11) 


CH3 


-54.2 (0.4) 


1.8 (0.6) 


258 (17) 


1.7 (0.2) 


97.0 (15) 


-62.8 (0.5) 


1.5 (0.9) 


927 (161) 


0.6 (0.1) 


97.3 (18) 


NH3 


-57.3 (0.3) 


0.9 (1.0) 


278 (31) 


5.2 (1.5) 


8.9 (1.2) 


-60.1 (0.2) 


0.4 (0.1) 


248 (32) 


2.8 (0.8) 


15.8 (3.6) 


CH4 


-56.3 (0.3) 


0.2 (0.04) 


674 (360) 


20 (8) 


5.8 (1.9) 


-60.1 (0.3) 


0.6 (0.3) 


668 (250) 


16 (13) 


99 (13) 


CS 


-55.2 (0.4) 


1.0 (0.3) 


224 (173) 


0.02 (0.03) 


100 (8) 


-59.4 (0.5) 


0.7 (0.3) 


249 (100) 


7.0 (3.0) 


5.5 (0.3) 



Note. 



The column densities listed are Ntot x C. 



molecules. 

Note that our quoted column densities are averages 
over the lines of sight to the continuum emitting mate- 
rial, or column densities in the absorbing gas multiplied 
by the covering factors. There are several possible inter- 
pretations for the covering factor parameters: the contin- 
uum source could in fact be partially covered, foreground 
or surrounding emission could veil the spectrum, or re- 
emission in the lines could fill in the absorption. In all of 
these cases, the column densities through the absorbing 
material are given by N/C, rather than N, the average 
column density over the line of sight in the partial cover- 
ing interpretation. Consequently, if we could determine 
the H2 column density along the lines of sight on which 
we see absorption we should compare that to our N/C 
values. 

4.1. H2 and CO From Previous Measurements 

Before discussing the results of our measurements of 
individual molecules, we attempt to derive H2 and CO 
column densities from previous measurements. The H2 
column toward NGC 7538 IRS 1 has not been directly 
measured, but it may be estim ated from the extinction 
measurements of lWillneil (|1976[ ) with the assumption of a 
normal interstellar gas-to-dust ratio. This method gives 
Nh2 = 6 — 9 X 10^^ cm~^. If the continuum source is 
non-uniformly covered, this column density represents 
an average over the source. Hence, it should probably 
be compared to the values we derive for average column 
densities, N, of the molecules we observe. However, the 
silicate optical depth includes extinction from all mate- 
rial along the line of sight to NGC 7538, most of which 
is unlikely to contain our molecules. This suggests that 
6 — 9 X 10^^ cm~^ is an overestimate of Nhs in the rele- 
vant gas. On the other hand, we suggest below that we 
may be observing absorption by gas in the atmosphere 
of a protoplanetary disk, and dust in such a region could 
be depleted by settling toward the disk midplane. This 
effect would result in an underestimate of Nh2 from the 
extinction. We give the ratios of our derived column 
densities to an H2 column density of 7.5 x 10^^ cm^^ in 
Table m 

It may be preferable to compare our column den- 
sities to the column density of CO, which has been 
observed in absorpt ion toward NGC 7538 IRS 1 by 
iMitchell et al.l ( 1990( 1. They obtain a column density of 
1.5 ± 0.3 X 10^^ cm-2 for ^^CO in cold gas at 25 K and 
1.4 ± 0.1 X 10^^ cm^^ in warm gas at 176 K. We would 



TABLE 4 

Abundances with respect to CO and H2 



Molecule 


N(X)/N(CO)* 


N(X)/N(H2)* 




10"=^ 


10-^ 


C2H2 


5.7 


7.7 


HCN 


6.9 


9.2 


HNCO 


0.5 


0.6 


CH3 


2.3 


3.1 


NH3 


8.1 


10.8 


CH4 


36 


48 


CS 


7.0 


9.3 



Note. — *N(CO) 
N(H2) = 7.5x10^2 



= 1.0x10^^ cm-2 and 
See text for discus- 
sion of whether these molecules are measured 
on the same lines of sight as the ones we ob- 
serve. We use A'^ listed in Table |3] 

not be sensitive to their cold gas. Presumably their warm 
gas includes both of our components, although our tem- 
peratures are generally higher than theirs. If we assume 
i^CO/^CO = 72, using the ^^C/^^C ratio we derived 
from our C2H2 observations, we derive the abundances of 
our molecules r elative to CO listed in TablelH (Note that 
IMitchell et all use i^c/iSC = 89.) However, we should 
note that we have measured one spectral setting includ- 
ing CO lines, while searching for OCS absorption near 
5 fim. The noise is rather high, but it is apparent that 
the ^'^CO lines are not well describ ed by moderate ly sat- 
urated Gaussians, as assumed by IMitchell et al.l . The 
^■^CO lines are deeper and broader than the lines of any 
of our other molecules, and have prominent blue-shifted 
shoulders probably tracing an outflow. We suspect that 
the lines are m ore saturated, and the CO column den- 
sity larger than IMitchell et al.l concluded. The different 
shapes and greater depths of the CO lines than our C2H2 
lines also indicates that CO probes different gas. Likely, 
C2H2 and the other molecules we observe in the mid- 
infrared are found only in unusual regions, whereas CO 
is distributed along the entire line of sight. 

In general, we conclude that the ratios of column densi- 
ties to those of either H2 or CO given in Table[4]are quite 
uncertain because the different column density tracers 
may be sensitive to different components along the line 
of sight. The column densities of the different molecules 
we observe should be much more comparable. 

4.2. C2H2 and HCN 
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Fig. 6.— The spectrum of NGC 7538 IRS 1 as seen by ISO 
is shown in black. A fit to the C2H2 and HCN fundamental Q- 
branches as well as a hot band Q-branch of C2H2 is shown in 
grey. The parameters are set from the TEXES high-resolution 
observations. The models gives a good fit to the ISO spectrum. 
Color version is available online. 



Since our initial search for molecules began with C2II2 
and HCN, the first models only included those two 
molecules. In these models, we derived similar tem- 
peratures and column densities a s derive d from ISO ob- 
serva tions (|Lahuis fc van Dishoe ck 2000; Boonman et al.l 
l2003f) . Once we include absorption from other molecules, 
HNCO in particular, the temperatures for HCN and 
especially for C2H2 are lowered significantly. Some 
HNCO lines overlap high- J lines of C2H2, which means 
that more hot C2H2 is necessary to produce the ob- 
served absorption when HNCO is not included in the 
model. In our final model including all of the ob- 
served molecules, our derived column densities of C2H2, 
TV = 5.8 X IQis cm-2, and HCN, = 6.9 x IQi^ 
cm~^, (see Table |3l) a re larger than those derived in 
iLahuis fc van Dishoeckl . A^C2H2 = 8 x 10^^ cm~^ and 
A'^HCN = 1.0 X 10^^ cm~^. Part of the difference is due 
to the differen t Doppler parameter used: 6 = 5 km 
is used in La huis fc van Dishoeckl . whereas, we derive 
6^1 km s^^. We derive lower temperatures, TC2H2 = 
190 - 230 K an d Trcn = 250 - 450 K than derived by 
ILahuis fc van D ishoeck. TcnUn = 800 K and Thcn = 600 
K, and lBoonman et al] (|2003l ). Tc^Hs = 500 K. We also 
find that we need a small covering factor to account for 
the saturated yet shallow lines. The level of saturation 
is difficult to determine from the ISO data. In order to 
verify that our fit to the TEXES data is also consistent 
with the ISO observations, we take our model at the ISO- 
SWS resolution and compare to the data (Figure |6|). The 
model matches the fundamental Q-branches for the two 
molecules as well as a hot band Q-branch of C2H2. 

C2H2 is expected to be in rotational local thermal equi- 
librium (LTE) since it has no permanent dipole moment. 
The fact that the two excited vibrational states of C2H2 
from which absorption was observed, and 1^5, are found 
to be populated near LTE at the C2H2 rotational temper- 
ature, or even slightly higher, provides information about 
the physical conditions of the gas. For both bands, the 
hot bands indicate they are slightly out of LTE since the 
vibrational temperatures (Ty) are higher than the rota- 
tional temperatures (Tr). For the + 1^5 - 1^4 bands, 
the Ty temperatures for the two components were 346 K 



and 271 K while for the 21/5 - 1/5 bands, they were 277 K 
and 227 K. The rotational temperatures for the V5 fun- 
damental are 225 K and 191 K for the two components. 
In LTE, the vibrational temperatures for the hot bands 
would equal the rotational temperature measured by the 
fundamental band lines. 

The 1^4 level has no allowed radiative transitions to the 
ground vibrational level, and none of the rotational levels 
of the ground vibrational level are coupled radiatively, in 
both cases due to the symmetry of the C2H2 molecule. 
This suggests that the populations of the ground and 1^4 
states should be set by collisions, and so should be in LTE 
at the same temperature. However, the 1^4 level can be 
radiatively pumped by absorption in the 7.6 /.tm 1^4 + vc, 
band followed by emission of photons in the 1^4 + 1^5 - 1^4 
band. The state must be populated predominantly 
radiatively for any plausible density, since the critical 
density is > 10"'^° cm~'^. From these considerations, we 
conclude that the brightness temperature of the radia- 
tion field seen by the observed gas must be close to the 
kinetic temperature of the gas. In addition, the fact that 
the measured rotational temperature for HCN is simi- 
lar to that of C2H2 indicates that the rotational levels of 
HCN are kept populated either by collisions, which would 
require a gas density > 10*" cm~ ^, or by infrared pump- 
ing th rough the 1^2 transitions. ILahuis fc van Dishoeckl 
()2000( ) conclude that either excitation mode is possible. 
The presence of C2H2 2i^5 - absorption requires in- 
frared pumping, although V4 + 1^5 - V4 absorption does 
not. 

We can attempt to estimate the extent to which emis- 
sion in our lines fills in the absorption, at least for C2H2, 
by using the 2v^ - lines as probes of the popula- 
tion. Our spectra require a j/5 vibrational temperature 
~300 K. We need to compare this number to the black- 
body brightness temperature of the background contin- 
uum rad iation. From the shape of the 8-13 /im spectrum, 
IWillneii derived a dust temperature of 330 K assuming 
the emitting dust is optically thick, or 370 K assuming 
optically thin silicate emission. He preferred the latter 
model, which gives a brightness temperature only when 
combined with an assumption about the solid angle sub- 
tended by the source. In either case, our vibrational 
temperature is only ~10% less than the background tem- 
perature, so that the source function is '^30% below the 
background intensity, and lines would be expected to sat- 
urate at '^70% of the continuum. However, the contin- 
uum spectrum from which .Willner derived a dust tem- 
perature is the average over the continuum source. If 
the gas we observe covers only a small fraction of the 
source, the brightness temperature of the radiation pass- 
ing through that gas could be much larger than the av- 
erage. 

Given the many assumptions and uncertainties in this 
calculation, it is difficult to make a strong statement. 
The level at which our C2H2 lines saturate, ~85% of the 
continuum, which we model with a small covering factor, 
could instead be due partly or entirely to emission by the 
absorbing gas. It is possible for the lines to saturate at 
85% of the continuum even with small covering factors 
(e.g., 24% and 6% for C2H2) since the line widths are 
narrower than our spectral resolution. Presumably the 
same statement applies to the other observed molecules. 
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4.3. NH^ 

From our infrared observations, we find that NH3 has 
two components: one at -57.3 km s^^ and the other at 
-60.1 km s""'^. Most of the NH3 is found in the first com- 
ponent, for which the temperature is ^ 278 K and the 
cohimn density is 5.2 xlO^^ cm~'^. The second compo- 
nent has a lower temperature (~248 K) and the column 
density is 2.8 xlO^^ cm-2(see Table E]). NH3 absorp- 
tio n toward IRS 1 ha s als o been studied in the radio 
by I Wilson et al.l (|1983[) and lHenkel. Wilson, fc JohnstonI 
(1198|)- They found A^NHa = 2 x 10^*^ cm'^ in gas with T 
= 170-220 K and Vlsr = -60 km s^^. Although their 
temperature and Doppler shift is in reasonable agreement 
with our bluer absorption component, their column den- 
sity is nearly 20 times ours. Without constraining the 
i5]\j^i4|sj isotopic ratio, our data allow an NH3 column 
as large as that from radio observations, but with this 
NH3 column -'^^NH3 lines should have been apparent in 
our spectra. It is possible that we sample only a fraction 
of the column observed at cm wavelengths as a result of 
larger dust opacity at mid-infrared wavelengths prevent- 
ing us from observing lines of sight through the densest 
regi ons of a knotty gas distribution. It may be notable 
that IWilson et al.l found, as did we, that although rela- 
tive depths of lines (in their case hyperfine components) 
require large optical depths, line shapes do not appear 
flat-bottomed as would be expected if they are optically 
thick. They suggest that the absorbing gas may be in 
small spectrally and spatially unresolved knots, each of 
which is optically thick, and which combine to make a 
line that does not have the shape of a thick line. Our use 
of two components in our fit may be a way of approxi- 
mating this situation. 

4.4. HNCO 

The velocities for the two components are -57.2 and 
-60.2 km s~^, with corresponding temperatures of 319 
and 171 K and column densities of 4x10^^ and 1x10^^ 
cm~^. Since the HNCO temperatures and velocities de- 
rived from our observations agree with the values for the 
other molecules, we can assume that HNCO is in the 
same gas. It appears th e region probed by HN CO is 
very small in angular size. IZinchenko et al.l (|2000f) detect 
HNCO emission in the radio and derive a beam-averaged 
column density of 1.1x10^** cm^^. This value is an order 
of magnitude lower than our observed column density. 
However, beam dilution in the radio observations can 
explain the difference in the column densities. 

4.5. CH3 

CH3 is detected for the first time toward warm, dense 
gas. The parameters derived for this molecule do not 
agree as well with the other molecules, although the 
weakness of the observed lines makes the derived param- 
eters rather uncertain. The separation between the cen- 
troid velocities for the two components, Vlsr — -54.2 and 
-62.8, is nearly twice the separation found for the other 
molecules. Also, the covering factor for both of the com- 
ponents is consistent with 100%, whereas as it is closer 
to 10% for the other molecules. In addition, the temper- 
ature of the second component is much hotter (>900 K) 
than the temperatures derived for the other molecules. 
We note though that the fact that the observed lines are 



of high excitation makes them particularly sensitive to 
a hot gas component, which may not have been noticed 
for other molecules. And the weakness of the lines makes 
the derived covering factor (which is inferred from rela- 
tive line depths indicating saturation) very uncertain. 

4.6. CHi and CS 

CH4 gas-phase absorption had be en observed to- 
ward the neighborin g protostar IRS 9 (|Lacv et al.lll991l : 
iBoogert et al.l [2003 l but had not previously been seen 
toward IRS 1. We find centroid velocities of -56.3 and 
-60.1 km s^^, which are similar to the other molecules. 
The temperatures are higher than for most molecules (^ 
670 K for both components) but have large uncertainties 
(see Table ^ . The best fit seems to indicate that one 
of the components covers the entire source whereas the 
other covers about 6% of the source. Of all of the ob- 
served molecules, CH4 has the highest column density. 
The implications of the high CH4 column density the 
chemistry are discussed in ij5.1l 

We detect only 6 lines of CS (see Table [H - full table 
available online only) . The lines are rather weak and pa- 
rameters are poorly constrained. Despite uncertainties, 
the parameters agree with those for other molecules. The 
centroid velocities are at -55.2 and -59.4 km s~^. The 
temperatures, Ti = 224 K and T2 = 249 K, are simi- 
lar to those found other molecules (except for CH4 and 
CH3). From radio studies of the envelope around IRS 1 , 
a CS abundance of ~ 10~^° is found (.Plume et al.llT997l : 
iMueller et all l2002t IShirlev et ani2003D. Ho wever also 
using radio observations, Ivan der Tak et all (|200Gf ) de- 
rive a CS abundance of ~ 10~*. Our abundance is about 
2 orders of magnitude larger than the largest estimate of 
the abundance in the envelope indicating that our obser- 
vations are probing material closer to the protostar. 

5. MODELS 
5.1. Chemical Models 

We now consider the implications of the observations 
for the chemistry of the observed material. Hot core 
models can be used to represent two of the scenarios 
discussed in !JTJ 1) material in a circumstellar disk and 
2) photo-evaporating knots of neutral molecular mate- 
rial. Hot cores are small (r < 0.1 pc), dense (riHa > 10'' 
cm~^) and hot (T > 100 K) regions associated with high- 
mass young stellar objects (jKurtz et al.ll2000( ). Chemi- 
cally, hot cores are identified by an enhanced abundance 
of fully hydrogenated molecules such as H2O and NH3, 
which are usually observed through rotational transitions 
at sub-millimeter and millimeter wavelengths. These 
abundances are enhanced with respect to the cold molec- 
ular envelope. A hot core is thought to form as a result 
of the protostar heating nearby material , which triggers 
ice su blimation from the grain mantles (ICharnlev et al.l 
[199^ . HNCO is another hot core molecule which has 
been found in the Ori on hot core with an abundance 
~ 10"^ relative to H2 dZinchenko et al.ll2000( l. 

Gas-phase HNCO is possibly coming from evapora- 
tion of grain mantles. The 4.62 /im feature observed 
toward some protostars has bee n attributed to OCN~ 
frozen in grain mantles ( e.g., iPendleton et al.l 119991 : 
Ivan Broekhuizen et al.112004 ). Recent ISO results of solid 
state features toward protostars show an OCN~ upper 
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limit of lO^'^ toward IRS 1 (|Gibb et al.ll2"00l . The 

total HNCO column derived in this work is 5.4x10^^ 
cm"^, a factor of 2 lower than the solid OCN" limit. In 
general, IRS 1 does not show many ice features compare d 
to the colder nearby source IRS 9 fsee lGibb et an[20nl . 
which has a column density N{OG'N~) = 1.2x10"^^ cm"^. 
If solid OCN" sublimates and captures a proton to form 
HNCO, the expected column density for HNCO would 
be comparable to the observed OCN^ ice toward cold 
lines of sight (like IRS 9). However, the gas-phase col- 
umn observed is at least an order of magnitude smaller 
than that of the solid state ion toward sources like IRS 
9. This indicates that when OCN" sublimates, a frac- 
tion of the ices goes into forming HNCO, while most of 
HNCO is destroyed in the warm gas-phase chemistry on 
short time scales. Further studies of how OCN^ reacts 
in the gas-phase are needed in order to understand the 
difference between ice and gas phase column densities. 

Similarly, CH4 is probably also evaporating from the 
grain mantles, and thus the gas-phase abundance is 
highly enhanced compared with the cold molecular enve- 
lope. Studies of the solid CH4 content toward IRS 1 show 
that the column density is 14 times less t han our ob- 
served gas-phase column ()Gibb et al.ll2004D . suggesting 
that the solid CH4 on grains has sublimated. In com- 
parison, IRS 9 has a higher content of solid CH4 tha n 
gas-phase CH4 (jBoogert et al.ll2004 iGibb et al.ll2004f ). 
If we consider th at IRS 1 is likely a more evolved ob- 
ject than IRS 9 (jElmegreen k Ladal Il977f l. the higher 
gas-phase CH4 toward IRS 1 also points to grain man- 
tle evaporation. The high abundance of CH4 means that 
daughter products such as CH3 and CH2 should be abun- 
dant. While the photo-dissociation of CH4 preferentially 
forms CH2 rather than CH3, our observations of CH3 
agree with this branching ratio. We predict that CH2 
is also present toward IRS 1 with a higher abundance 
than CH3. It is uncertain where C2H2 is formed in the 
gas-phase chemistry or is evaporating from grain man- 
tles. However, most models do not predict the observed 
column densities from formation in gas-phase chemistry 
only. 

From our observations we see a high abundance of the 
parent molecules NH3 and CH4. Hi gh abundances of 
the ot her molecules are also observed. iNomura &: Millail 
()2004|) present a time dependent chemical evolution of 
molecules in a hot core. The models begin at the time 
when the protostar turns on = yrs). Table [5] shows 
the abundances with respect to CO at an age of 10^ 
years. The hot core model seems to underpredict abun- 
dances with respect to CO of all our observed molecules. 
However, as mentioned in the CO column den- 



TABLE 5 

A''(X)/A''(CO) FROM CHEMICAL MODELS 



sity found bv iMitchell et al] ( 1990f) may not be the ap- 
propriate number with which to compare our observa- 
tions. If the CO column along the lines of sight con- 
taining o ur molecules is larger than the column mea- 
sured by IMitchell et al.l . then perhaps the abundances 
from the model will more closely resemble the observa- 
tions. We can also compare the relative abundances be- 
tween our observed molecules, avoiding the uncertainties 
in whether the CO and dust absorption trace the lines 
of sight that our molecules are found in. The observed 
and predicted values for A^(CH3)/A^(CH4) are similar, 
~0.01. For the other molecules the ratios from the mod- 
els at an age of 10'' years are very different from the 



Molecule 


Observed 


Hot Core Model" 


Disk Models' 




(10-3) 


(10-3) 


(10-3) 


C2H2 


5.7 


0.88 


0.06 - 0.1 


HCN 


6.9 


0.02 


16 - 24 


CH3 


2.3 


0.016 


0.02 - 1.6 


NH3 


8.1 


0.76 


1-6 


CH4 


36 


1.6 


7-35 


CS* 


7.0 


0.0031 


0.03 - 0.07 



References. — " INomura fc Millarl2004l and H. Nomura, 
private communication 2004 
INguven et al.ll200^ 
Note. — The model values were taken at 10* years. 

observed values. The values for the models assume that 
we can see all the material and comparing the model re- 
sults directly with the observations without accounting 
for line-of-sight does not make for a good comparison. 
Figure [7] shows the column density ev olution for 6 of the 
detect ed molecules as determined by iNomura fc Millail 
(|2004l ) corrected for the observed beam (priv. comm.). 
For each molecule, the observed range of column densi- 
ties is shown in the horizontal shaded area. This hot core 
model can predict the observed column densities for 4 of 
the 6 molecules plotted. However, the ages indicated by 
each molecule do not give a consistent overall age. For 
CH4 and HCN the model overpredicts the observed col- 
umn densities. The other molecules suggest ages between 
2x lO'^ and 2x 10^ years. While ages vary by three orders 
of magnitude, it should be noted that the model has not 
been optimized to be an accurate representation of this 
source. With chemical models that more closely repre- 
sent the source, it may be possible to determine a more 
consistent age. 

If we consider the disk scenario, we can also com- 
pare our observat ions to a chemical m odel of a photo- 
evaporating disk. iNguven et al.l ()2002[ ) present a chem- 
ical model for a disk around a 10 M0 star. The model 
consists of a flared disk with a photo-evaporating layer 
on the surface. Table[5]compares our derived abundances 
with respect to CO for the various molecules to the pre- 
dicted values from disk models. The summarized results 
of models taking the midplane temperature for the tem- 
perature of the layers below the photo-dissociation re- 
gion (PDR) are discussed here. The difference among 
the three models is the variation in the X-ray ionization 
rate, C (IC = 1.3x lO^i^s-i), from 1 to lO^C- Tabled 
shows that the abundance for C2H2 is underpredicted. 
In the chemical model, C2H2 is produced only by gas 
phase reactions. One possible explanation is that C2H2 
is formed as ice on grain mantles. The abundance of 
NH3 is also underpredicted. The model with the highest 
X-ray ionization rate produces CH3 abundances close to 
the observed value. For that model, 7V(CH3)/iV(CH4) 
is 4 times larger than the observed ratio. These models 
do not fit our observations very well but can serve as a 
guide to t he expected ab undances in such a scenario. The 
model by INguven et al.l has not been modified to try to 
fit data. So, adapting the model for this object may give 
predictions closer to the observations. The star in NGC 
7538 IRS 1 is believed to be a 30 Mq star. The different 
radiation field may affect the resulting chemistry. 
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Fig. 7. — These plots illustrate the determination of ages 
based on the column densities of the various molecules compared 
with time dependent ch emical model. The model data are from 
INomura &: Millar! II2004I ) . The horizontal lines indicate the range 
in column densities derived in this paper. The vertical lines indi- 
cate the age range corresponding to those column densities. Note 
that for CH4 and HCN the model overpredicts the column density 
at all ages. 

iDotv et all (|2002D present a model of chemical evolu- 
tion of the envelopes of massive protostars using AFGL 
2591 as an example. From that model, they conclude 
that enhancements of C2H2, HCN and CH4 are possi- 
ble at late times in the evolution (> 10*^ years) and at 
high temperature (T ~ 800 K). Our derived tempera- 
tures indicate cooler material, yet the molecules seem 
to be enhanced in the line of sight. Alternatively, these 
species could be produced at high temperature in the 
inner disk, and subsequently brought outward by the 
disk wind. Including UV and X-ray radiation in the 
model affects t he ch e mistr y especially relating to HCN. 
iStauber et aH (l2003 . [200I find that' increasing the UV 
and X-ray flux can lead to enhanced HCN starting by the 
ionization of N2. Thc_column densities for HCN and CS 
from the fStauber et al.i model of protostar AFGL 2591 
agree with our observed abundances while their column 
densities for C2H2, CH4, and NH3 are lower than our 
observed values. 

Hot core and disk chemistry models predict the en- 
hancement of molecules such as C2H2, CH4, and NH3. 
However, the observed abundance of C2H2 is higher than 
what models predict from warm gas-phase chemistry (see 
Table [5]). This indicates that C2H2 is probably frozen 
on du st grains and sublim ates along with molecules like 
CH4. iBoudin et al.l ()1998D studied the solid features of 
C2H2, especially when mixed with H2O or CO. They 
found that the C2H2 features broaden substantially when 
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Fig. 8. — A cartoon representation of the IRS 1 region. The 
absorbing material containing the molecules observed can either 
be located in the disk or in the hot knots, while the continuum is 
coming from the photo-evaporating layers of the the inner rim and 
surface of the disk. Flaring is not to scale. The continuum can 
also be from other parts of the disk for which lines of sight do not 
cross the material with our molecules. There may also be more 
extended dust emission from the envelope that contributes to the 
continuum but does not contain our molecules in high abundances. 
In this cartoon, North roughly points toward the top of the page. 

mixed with H2O and to a lesser extent when mixed with 
CO. They compare the laboratory data to ISO data for 
the colder neighbor, IRS 9, which resulted in an upper 
limit of 8 X 10^^ cm~^ for the column density of solid 
C2H2. This upper limit is consistent with the observed 
gas-phase column density seen toward IRS 1. So, it is 
possible that solid C2H2 is sublimating from grain man- 
tles as protostars heat the environment. The observed 
CH3/CH4 ratio agrees with the branching ratio for the 
destruction of CH4. Based on these results CH2 should 
also be very abundant. 

5.2. Physical Models 

We now attempt to construct a physical and geomet- 
rical model of NGC 7538 IRS 1. Within the possi- 
ble scenarios prese nted from va.rious r adio and infrared 
observations (e.g.. iMinier et all 120011: iLugo et all 120041 : 
iDe Buizer fc MinieHl2005tlKraus et "311120061 ) . we propose 
a scenario in which the molecular absorption presented 
here comes from a circumstellar disk. We will examine 
other possibilities first. Figure [8] depicts the possible sce- 
narios allowed by the available observations. 

We can quickly rule out a simple model in which the 
absorbing molecular gas is in the foreground molecu- 
lar cloud and not closely associated with IRS 1. Our 
observations require the gas to be much hotter (T ~ 
300 K) and denser (riHa ~ 10^ cm""^ to maintain rota- 
tional LTE of HCN out to J — 21) than is found away 
fr om luminous sourc e s in m olecular clouds. According 
to Ivan der Tak et ahl ()2000r ) , temperatures do not reach 
the observed values until you get to within 400 AU of 
the central star. In addition, the infrared radiation field 
must have a brightness temperature ~300 K to populate 
the C2H2 I's level sufficiently to account for the observed 
2i^5-i^5 absorption. Undoubtedly, the absorbing molec- 
ular gas is in close proximity to the infrared continuum 
source IRS 1. 

It is not quite so easy to rule out a model in which 
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the absorbing molecules are in boundary region between 
the envelope around the IRS 1 hypercompact H II re- 
gion and the outflow. This gas could be compressed 
by the ionized wind, and if it is as close as 0.1", or 
280 AU, from a 10^ L© source it would have a tempera- 
tu re near 300 K. The tempe rature structure determined 
bv Ivan der Tak et al.l (|2000f) indicates that indeed tem- 
peratures range from 200 - 400 K at radii of 220-380 
AU. However, interaction with the ionized wind, which 
has a velocity ~100 km s~^, would be expected to accel- 
erate the gas, causing br oad, blue-shifted absorption (see 
Ivan der Tak et al1l2000( ). In contrast, the observed hues 
have widths of < 8 km and have centroids within a 
few km s^^ of those seen in surrounding molecular gas. 
Another argument against the presence of the observed 
gas in an envelope around IRS 1 is the fact that net ab- 
sorption is seen. Since the vibrational temperature is 
comparable to the brightness temperature of the contin- 
uum radiation, emission lines would be seen if the molec- 
ular gas had a larger extent than the continuum source. 
If the lines arose in molecular shell surrounding the hy- 
percompact H II regi on, this probably would be the case. 

On the other hand, ICampbeOI ()1984[ ) proposes that the 
centimeter continuum from IRS 1 results from partially 
ionized material in an outflow. The centimeter con- 
tin uum emission has been spatially resolved into knots 
bv iGaume et al.l ((T995i ). who suggest that the emission 
comes from photo-evaporation of knots of neutral molec- 
ular material (see F ig. [8]). The CII3OII masers seen by 
iMinier et al.l (|2000( l also seem to trace the knotty struc- 
ture (in addition to the disk described in [JIJ probed by 
the cm continuum emission. The knots may be mate- 
rial stripped from the disk. The stellar wind would then 
blow on these knots. However, if our molecules are in 
these knots we would expect large velocity difference be- 
tween the knots resulting in broad lines, which are not 
observed. If we compare the velocities we derive for the 
two components of our molecules to those observed for 
the CH3OH masers, we find that they match the two 
velocities found for the A cluster of masers. This indi- 
cates that the infrared observations may be probing the 
same gas as t he CH3OH masers. This scena rio is further 
supported bv lKraus et al.l ()2006l ). However. iKraus et al.l 
also find that the interpretation of the masers tracing a 
disk also agrees with their observations. 

We propose that the absorbing material is in 
a circumstellar disk. There is evidence for a 
disk close to edge-on from other observation s and 
models (e.g., [ Pcstalozziet_alJ 120041: iLugo et al.l 12004 
iDe Buizer fcMinieri,200'5[ r In this situation, hot dust in 
the surface layers or the inner rim of a flared disk could 
provide the continuum source. With a nearly edge-on 
orientation, radiation from the inner region could pass 
through the outer disk atmosphere, where absorption 
lines could be formed. Emission from some of the dust 
on the surface of a slightly inclined disk or in an outflow 
could be observed directly, accounting for the dilution 
of the spectrum that we modeled with a small covering 
factor. This scenar io is similar to the one described by 
iLahuis et all (|2006[ ) to describe C2H2 and HCN absorp- 
tion toward the low-mass protostar IRS 46 in Ophiuchus 
(YLW 16B). Detailed modeling of the disk structure of 
IRS 46 showed that it was possible to have the absorb- 



ing material in the disk illuminated by the continuum 
from the surface. Figure |8| is a cartoon of our proposed 
picture, where the orange rays represent the continuum. 
Some of the continuum radiation passes through the ab- 
sorbing outer disk. In this scenario, dust settling is likely 
to occur, allowing for observations of large columns of gas 
without too much dust. 

6. CONCLUSIONS 

Here we present a summary of our conclusions: 

1. NGC 7538 IRS 1 has a rich spectrum in the mid- 
infrared. At high spectral resolution, lines of previously 
observed molecules such as C2H2 and HCN are seen as 
well as weak lines of molecules such as HNCO, CH3, NH3, 
CH4, and CS. We have presented the first infrared detec- 
tion of interstellar HNCO. This is also the first detection 
of CH3 toward dense gas. 

2. The data show shallow yet saturated lines indicat- 
ing that either the absorbing gas does not fully cover the 
continuum source or there is emission filling in the ab- 
sorption. The observed ortho:para equivalent width ratio 
for C2H2 is ^ 1.5 instead of the expected value of 3. 

3. The rotational temperatures for the various 
molecules range between 200 and 400 K with the excep- 
tion of CH3, which has one hot component (T ~ 900 K). 
From the derived temperatures, we fi nd that the gas is 
withi n ~ 400 AU of the central star ()van der Tak et aTl 
l2000[) . 

4. We suggest that the material traces a close to edge- 
on circumstellar di sk as determined f rom other observa- 
tions fe.g.. iDe Buizer fc Min icr 2005). However, it can 
also be tracing hot, knotty material close in to the star 
(see Figure E]). 

5. Because it is possible to use hot core chemistry 
to describe chemistry in a disk, the different chemical 
models do not help identify the physical location of the 
absorbing gas. However, the models do seem to indicate 
that NGC 7538 IRS 1 is an evolved protostar with an 
age ^ 10^ years. 

6. Our observations help constrain the chemistry in 
massive protostars. The abundances for all the molecules 
are enhanced, except for possibly CS (see discussion in 
§ 14. 6p . Species known to be present in icy mantles such 
as CH4 and NH3 have high column densities suggesting 
that they have recently sublimated from grain mantles. 
The high abundance of the daughter molecules such CH3 
gives constraints on the gas-phase chemistry happening 
after sublimation of ices. HNCO also provides some con- 
straints about the presence of OCN~ on dust grains. It is 
unclear whether C2H2 is a parent molecule or a product 
of warm gas-phase chemistry. For molecules like HCN, 
UV/X-ray radiation is important in producing the ob- 
served abundances. 
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HNCO Selection Rules 
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APPENDIX 
SPECTROSCOPY OF CH3 AND HNCO 

CH3 is a nearly planar molecule and thus do not expect to see inversion transitions as observed in NH3 . We observe 
the out-of-plane vibrational mode (similar to the umbrella mode in NH3). We observe splitting due to spin-rotation 
interactions. For symmetric-top molecules, the splitting of a rotational level is given by 

Az. = (iV + l/2){e':; - (e':;' - tJ^K^/lNiN + 1)]}, (Al) 

where e\l and t \.\ are the spin-rotation coupling constants for the B and C axes, respectively, at a given vibrational 

state, v. We use Equation ([Alll in deriving the splitting for the lines we observed. The values for e'^;,' and e were 
taken from I Yamada et aTl ([1981,). 

HNCO is a quasilinear, nearly symmetric-top molecule. The moment of inertia about the figure axis is small and 
the moments of inertia about the other two axes are much larger (~ 100 times larger) and about equal to each other. 
Because it is a slightly asymmetric molecule, the levels with K 7^ are split into two components. iSteiner et al.l ()1979f ) 
refer to these levels as upper (U) and lower (L). Table [5] shows the selection rules for HNCO. 

In Figure m the HNCO lines labeled R(J)1U correspond to lines where A/f = and A J = 1 and K = 1 upper. 
Likewise, the lines labeled R(J)1L are lines where lAK = and A J — 1 and K = 1 lower. Because the splitting 
between the two levels is small the lines overlap in the spectral region shown in Figures [3] and [H Some lines from the 
upper and lower K states are blended together and give the appearance of double peaked lines. 
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